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ABSTRACT 

This thesis addresses the problem of blind adaptive equalization of finite 
impulse response (FIR) and infinite impulse response (IIR) channels (their 
characteristics are unknown) in the fractionally sampling scenario. An adaptive 
equalizer is used at the receiver to compensate the time dispersion induced by noisy 
communication channels and eliminate the effect of Inter-Symbol-Interference (ISI). 
In other words, the overall our system model, which is a cascade connection of the 
channel and equalizer, provides nearly an ideal transmission medium that the 
information source signals can be sent through. Due to this and rely only on 
probabilistic and statistical properties (Second Order Statistics (SOS) which has most 
communication channel information) of the received signals, the unknown input 
information signals can be recovered successfully. 

Various blind adaptive algorithms are discussed throughout this thesis. 
Simulation results are presented by evaluating the mean square symbol error (MSE) 
of these techniques to study their performance behavior in blind channel equalization 



in 



concept. These algorithms operate blindly in the practical situation, and they can 
achieve a complete adaptation without the aid of a training sequence, desired 
response, which is either impractical or very costly. The parameters of equalizer are 
updated in a recursive way with each single output measurement. Finally, the 
performance comparisons are realized to show which algorithm is more efficient and 
robustness to noisy channel model (three classes of channel model are used through 
thesis's simulation). The aim of this thesis is to improve the performance of a wireless 
communication channel using various blind adaptive equalization algorithms through 
computer simulations. 
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CHAPTER I 
INTRODUCTION 

1.1 Introduction and Background Information to Blind Adaptive Equalizer 

Communication systems have been growing during the past 40 years, and 
their applications have been used in many electronic products. In these days, radios, 
televisions, mobiles, and computer terminals with Internet comprise an essential 
element in our life. They all have ability of providing fast communications from 
every corner of the globe. We notice that there is an almost unlimited or endless 
amount of applications relying on the use of communication systems. 

Any digital communication system involves transmission of analog signals 
through channels. These channels are usually dispersive mediums that introduce 
delay or memory in the received signal, which in turn spreads the symbols over time. 
This spreading of symbols induces a distortion known as inter- symbol- interference 
(ISI) [11]. This kind of distortion is undesirable because it makes communication 
systems less reliable and in order to maintain our systems reliable performance, it 
must be removed at the receiver by equalization. In many high speed data rates 
bandlimited digital communication systems, it is shown that ISI appears in all 
communication channels such as twisted pairs cables, coaxial cables, fiber optics, 
satellite, microwave and radio channels. This interference between the adjacent 
symbols corrupts the detection of those symbols since their spreading can corrupt the 



time spacing. Increasing ISI causes to increase the error probability at the receiver. 
Therefore, under strong ISI, a simple memoryless decision device may not be able to 
recover the original data sequence [7]. 

ISI needs to be suppressed at the receiver to keep our systems in reliable 
performance. Equalizers are used as a solution at the receiver to compensate the 
distortion (ISI) occurring during transmission. Equalizers try to extract the 
transmitted symbol sequence by countering the effect of ISI. Therefore, they improve 
the probability of correct symbol detection [16]. 

One of the earliest channel equalization methods is achieved by sending a 
training sequence to estimate the linear filtering dispersion characteristics of the 
nonideal channel (channel identification). Note that the nonideal channel 
characteristics are usually not known priori, and the training sequence was already 
known at the receiver. The training sequence is required to be sent periodically since 
an adaptive channel equalizer relies on it, which occupies much of the bandwidth, 
resulting lower communication link efficiency. For this reason, a self-recovering 
approach that does not need a training sequence is invented. Such approach is termed 
'blind'. This idea of self recovering blind equalization was first introduced by Sato in 
1975 [7]. 

In recent few years, many different blind algorithms have been introduced, 
and it has been growing to the idea of the so-called 'blind problem' [4], [5], [14], [17] 
[28] and [29]. In general, blind processing can be defined as a digital processing of 



unknown signals that are sent through a linear channel without any knowledge of its 
characteristics and additive noise. In wireless communications systems, the blind 
signal processing plays an important role in combating frequency selective fading and 
inter-symbol interference (ISI) [12]. Blind adaptive equalization methods treat the 
problem of recovering the transmitted input sequence by using only the received 
output signal without knowledge of the channel characteristics. Since then, blind 
channel equalization has obtained great research interest because of the fact that 
receivers start equalization with no need of transmitter assistance (sending training 
signals). Consequently, a better channel bandwidth efficiency is obtained. 



1.2 Thesis Overview 

This thesis considers the problem of blind adaptive equalization in 
communication systems. It starts presenting some introduction to general 
communication systems and their applications in our daily life. The introduction also 
explains how the distortion (ISI) is generated and due to its effects on our system's 
performance, it has to be removed by using equalization. It introduces one of the 
earliest equalization algorithms (using training sequences) and its drawbacks with 
respect to valuable performance and throughput, causing to develop the so-called 
blind equalization techniques. Definition of blind processing and its advantages is 
also demonstrated in a simple way. Chapter 2 describes the channel model used 
throughout this thesis. Fractionally sampling scenario is considered in our channel 
model that is SIMO scheme. It also presents the complete system model (channel 
with linear equalizer) and their vector forms representation that will be useful for 
derivations the blind algorithms that are discussed in chapter 3. This chapter ends 
with discussing the conditions to achieve a perfect equalization. 

Chapter 3 begins with the definition of adaptive filter and its applications, 
including channel equalization. It introduces blind adaptive SIMO channel 
equalization method (using only second order statistics (SOS) of the received signals 
to design the equalizer) such as Direct Blind MMSE Equalizers. This chapter also 
presents the calculation methods for blind equalizer including derivations of 
Recursive Least Square (RLS), cyclic Least Mean Square (cyclic LMS), Least Mean 



Square (LMS), Linear Prediction (using GP technique) and Recursive Extended Least 
Square (RELS). 

Chapter 4 gives simulation results of the equalization techniques to examine 
their performance behaviors. One of the performance measures is Mean Square Error 
(MSE) that can be used to compare the equalization algorithms. The channel model 
classes used through simulations are also described in this chapter. The rest of the 
chapter is graphs, results, and comparisons between all of these algorithms with 
respect to their MSE using different classes of channel model. Chapter 5 gives a brief 
conclusion to the overall concepts included in this thesis. 



1.3 Mathematical Notations 

w c (t) Continuous-time input information signal 

w(n) Input discrete-time sequences 

w(ri) Equalizer output 

8(t) Dirac delta function 

h c (t) Continuous-time impulse response of composite channel 

h(ri) Discrete-time equivalent of h c (t) 

/i (n), /ii(n) Finite-length impulse responses of the two subchannels 

h Vector form of channel impulse response 

h , h-L Vector form of the two subchannels coefficients 

v c (t) Continuous-time additive noise 

v (n) Discrete-time equivalent of v c (t) 

x c (t) Continuous-time channel output 

y c (t) Continuous-time channel output with additive noise 

y(n) Received discrete data signal, input to equalizer 

y (n), y t (n) Received discrete data signals of the subchannels outputs 



g(n) Discrete-time impulse response of linear equalizer 

g Vector form of linear equalizer 

go, gi Vector forms of the two subequalizers 



T s Pulse duration 



* Convolution 

(. )* Complex conjugate 

(. ) T Transpose 

(. ) H Hermition or the conjugate transpose 

(. ) + Pseudoinverse 

H T Channel convolution matrix, Sylvester or Block Toeplitz 

matrix 

C2 W (n),C2 V (n),C 2y (ji).. Correlations of the stationary input, noise, and channel 
noisy output 

/ Cost function 

A 

V/ (JV) The instantaneous approximation to the gradient of the cost 

function J Q 

7(g) Prediction error variance 

II. II Norm or absolute value 



E[.] Expected value 

&w Input variance 

&% The channel noise variance 

X Forgetting factor 

\i, a The small step size 

y Constraint value 

q~ x Unit delay operator 

Biq' 1 ), C(q _1 ),i4(q _1 ).. Polynomials of IIR channels 

F(q _1 ) Common factor (minimum phase) 

y(i + 1) One step-ahead prediction of x 1 (i + 1) 

A 

y(i + 1) A posteriori prediction 

i?(q~ 1 ),5(q~ 1 ),D(q~ 1 ) .. Polynomials of IIR equalizer 



A 

6 (i) An estimate of unknown 9 



I Identity matrix 

p An arbitrary finite positive scalar 

MSE Mean square error 



P Roll-off factor 

r c (t,/3) Raised-cosine channel 

8 A parameter defining the size of the channel modeling 

error 

i.i.d Independent and identically distributed 

FIR, IIR Finite impulse response, Infinite impulse response 

SIMO Single input-multiple output 

SOS Second order statistics 

QAM Quadrature amplitude modulation 

LTI Linear time invariant 



CHAPTER II 

MODEL OF COMMUNICATION SYSTEM AND 
MATHEMATICAL FRAMEWORK 



2.1 Introduction 

A general and basic model for a communication system is comprised of digital 
signal that is transmitted by a transmitter through an analog channel to a receiver (see 
Figure 2.1). 



Source Info. 



Transmitter 



Channel 



Receiver 



Received Info. 



Figure 2. 1 General model of a communication system. 

Source information may be either a digital or an analog signal. The transmitter 
is fed from source information, which for the purpose of this thesis we assume is a 
binary. However, when the information source is an analog signal, transmitter is 
extended to have a so-called source encoder to convert the analog signal into a 
sequence of binary digits. This sequence is passed through a channel encoder to 
introduce a redundancy in the binary information sequence. This redundancy in the 
binary information is used at the receiver to overcome the effect of noise and 
interference in the transmission signal through the channel. The binary sequence at 
the output of the channel encoder is passed to a digital modulator which does the job 
of mapping the binary information sequence into waveforms to be sent through 
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communication channels. Figure 2.2 illustrates basic elements in a digital 
communication system model. The main digital modulation scheme used in the 
analysis and simulations throughout this thesis is a 16-QAM (quadrature amplitude 
modulation) since the symbols are usually complex with in-phase and quadrature 
components. Non ideal analog media such as telephone cables and radio channels 
typically distort the transmitted signal. During transmission, the essential features of 
the transmitted signal are corrupted in a random manner. 
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Figure 2.2 Basic elements in a digital communication system. 

At the receiver end of a digital communication system, the digital 
demodulator processes the channel corrupted signal waveform to convert the 
waveform into a sequence of number that represents the transmitted signal. Then, this 
sequence is passed through a channel decoder, which attempts to recover the original 
information sequence. How well the original signal is recovered depends on many 
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factors such as transmitted signal power, characteristics of the channel, the amount of 
noise and nature of interference. 

QAM 16 and higher is widely used in communication systems, and it 
combines both phase shift keying (PSK) and amplitude shift keying (ASK) 
modulations. A simple concept of input-output communication system using M-QAM 
is that the transmitter generates a sequence of complex-valued random input 
data (w(Z)}, each element of which belongs to a complex alphabet A (or 
constellation) of M-QAM symbols. The data sequence (w(/)} is sent through a 
baseband equivalent unknown complex linear time invariant (LTI) channel whose 
output x(t) is observed by the receiver. The function of the receiver is to estimate the 
original data {wQ)} from the received signal x(t) [7]. QAM 16 will be used in this 
thesis because its constellations allow for more bits per symbol and thus more bits per 
hertz. 

2.2 Channel Model 

In a digital communication system, the transmitter modem takes R bits of 
binary data at a time and encodes them into one of 2 R analog symbols for 
transmission (at the signaling rate) over an analog channel. At the receiver the analog 
signal is sampled and decoded into the required digital format. 

Assume that at the transmitter modem, the I th set of R binary digits is mapped 
into a pulse of duration T s seconds and an amplitude w(l). Thus the modulator output 
signal, which is the input to the communication channel, is given as 
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w c {t) = Ziw{l)S{t-lT s ). (2.1) 

where w c {t) is the continuous-time information signal, 8{t) is a pulse of duration T s 
with an amplitude w(Z). 

The unknown complex- valued LTI communication channel with the impulse 
response h(t) composited with the known transmit and receive filters is shown in 
Figure 2.3. The overall composite channel h c (t) is finite impulse response (FIR). 
Additive noise v c (t) that is assumed to be stationary and uncorrelated with the 
information symbols w(ri) can be included in this model. 

v c (t) 



Wc(t) 



h c (t) 



x c (t) 

► 






y c (t) A y( n ) 



p 



Figure 2.3 Fractionally sampled (P = 2 or more) communication system. 

The composite channel (whose impulse response is h c (t)) is non ideal, which 
means h c (t) is nonzero for t =£ 0. An integer P in Figure 2.3 denotes the amount of 
fractionally sampling (oversamp ling). The channel output can be modeled as the 
convolution of the input signal and composite channel response. The signal at the 
sampler is: 



ydt) = (w c (t) * h c {t)) + v c (f), 
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_ oo w c (r)/i c (t - t) dx + v c (t). (2.2) 
substituting (2.1) into (2.2), we get 

ydt) = Zr=-coL>(0«(T - Zr s )/i c (t - r)dr + v c (t). (2.3) 

y c (0 = Zi w(Z)/i c (t - Zr s ) + i; c (t). (2.4) 

If y c (t) is sampled at t = s /n, the received discrete data signal (y(ri) = 
y c (t)fort = n s /p)is 

y(n) = Ij w(/)/i c r r V P - W s j + v c ( nTs / p ). (2.5) 

Considering w(n) and y(n) are the input and output discrete-time sequences, 
respectively, it is convenient to rewrite (2.5) as an equivalent discrete-time system 

y{n) = ir=-oo w(/)/i(n - IP) + v(n), 

= x(n) + v(ri). (2.6) 

where h(n) and v(n) are the discrete-time equivalents of h c (t) and v c (t), 
respectively, x(n) is the discrete-time equivalent of the noise free received signal, 
and index n can be any integer (n = •■• , —3, —2, —1, 0, 1, 2, 3, •■• ) [11]. 
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In the single input-multiple output model, an integer P > 1 of measurements 
are performed for each transmitted symbol, provided the continuous-time observation 
is oversampled compared with the transmitted symbol rate, leading to create a 
fractionally spaced scenario. Throughout this thesis, it is assumed that P = 2 to make 
our model in the oversampling scheme (see Figure 2.4). 



w(n) 



h (n) 



hi(n) 



v (n) 



x o(n) 



r\ ^ y (n) 



Vi(n) 



*ifrO 



■e 



y t {n) 



Figure 2.4 SIMO multichannel model (P = 2). 



now, (2.6) can be rewritten in fractionally spaced system as follows 



yM) = E£-oo w(.l)h t (ri -0 + v t (ri) for i = 0, 1 



(2.7) 



with 



y (n) = w(n)* h (n) + v (n), 



Z^-coW(Z)/i (n-Z)+Vo(n). 



(2.8) 



and 
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yiOO = w(n) * /ii(n) + v-lOi), 

= Sr=-oo w(0/ii(n - + ^(n). (2.9) 

Note that /i (ft) and /^(n) are finite- length channel impulse responses of the 
first and second channels, respectively, at rate -y^ . 

2.3 Vector Representation 

A vector representation is very important and useful in our system analysis 
and equalizer implementation because it makes our finite-length subchannels 
/i (n) and h-^fri) easy to understand and modify. The input, noise and output such as 
w(n),v(n) andy(n) can be also expressed in vector form. The index n in these 
signals is usually infinite length. However, it will be finite-length when it is used as a 
vector form. Assume that the finite-length channel h(n) is causal with length 2L h 
such that each subchannel has a length L h . Considering the following vectors: 

y(ri) = [y {n) y^n)] 7 . (2.10) 

v(n) = [^o(ri) v 1 (ii)] T . (2.11) 



h = 



7i (0) /i (l) - /ioO-fc-1) 
Ai(0) /i!(l) ... /iiCLfc-1) 



(2.12) 



with 



h o = [/i o (0) /i (l) - /lo(^-l)]. (2.13) 
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and 



then 



h t = [^(0) h t (l) ... h t (L h - 1)]. 



y(n) 



y (ri) 
Viin) 



7i (0) /i (l) - 
AiCO) /lid) ... 

w(n) 
w(n — 1) 

w(n - (L h - 1)) 






+ 



i7 (n) 
v^ri) 



(2.14) 



(2.15) 



or 



r^h-i 



y(n) = 



Vo(n) 
yi(n) 



} w(n - l)ho(l) + v (n) 

1=0 
L h -1 

y w(n — /)/ii(/) + ^(n) 



i=0 



'^h-i 



= Llt w(n-lMl) + v(n). 



(2.16) 



where h is a vector form of the two subchannels, h and hj represent the vector form 
coefficients of the first and second channels, respectively. The noisy output vector of 
1-input/ 2-output causal 2L h FIR channel filters driven by the input sequence is y(n). 
Note that (L h — 1) is the subchannel polynomial order. 
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2.4 Overall System Model with Equalizer 

A finite impulse response filter g(ri) as an equalizer with vector parameters g 
is applied to our channel model to make our system into a complete system form (see 
Figure 2.5). The goal of the equalizer is to remove the distortion caused by channel 
ISI. 



v (n) 



w(ri) 







x (n) ^ 


V°( n ) . 








h (n) 


0oOO 






*\u 










v x {n) 








x t (n) ^ 


Syi( n ). 






hi(n) 


0i 00 




1 ► 


V 


V 





w(ri) 



Figure 2.5 Multichannel system models with equalizer g(j\). 

Assume that the finite-length equalizer vector g can be defined in a same 
manner as h (as mentioned in previous section) with length 2L g such that each 
subequalizer has a length L g . The following vectors are defined 

g = [0o (0) flr (l) 0o(^ - 1) 

0i (0) 0i(l) 0i(^-l)] T . (2.17) 

with 



go = [0o(O) 0o(D 0o(^ - l)f. 



(2.18) 
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and 



gi = [0i(O) 0i(D ^i(^-i)] ] 



(2.19) 



where g and g! represent parameter vectors of the first and second equalizers, 
respectively. Let w(ri) be the output signal of equalizer, and then the input-output 
relationship can be written in the following form 

w(n) = w(ri)* [h (n) * g (ri)] + v (n) 



+w(ri) * [/ii(n) * 0i (n)] + Vi(n), 



(2.20) 



w(n) = w(n) * [ 2^ g (l)h (n - /)] + v (n) 



i=0 



+w(n) * [Z^fo 0i(O/ii(n - 0] + ^i(n). 



(2.21) 



To simplify (2.21), let's consider our system in the absence of noise and define 
H T that is associated with vector h to be a channel convolution matrix with the length 
(l h + L g - l) X PL g as below 



W = 



/i (0) 

fci(0) h (0) 

/i (l) ^(0) 

/iiO-fc-l) /i (L ft -l) 

6 
o 



/i (0) 








Jllo* - 1) 



(2.22) 
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The convolution matrix H T is also called Sylvester or block Toeplitz matrix 
[15] of the two subchannels, and P is equivalent to linear filters that the input signal 
w(n) passes through (P = 2). Alternative representation of (2.22) is found in some 
books and is written in the following form 



H T = 



ho(0) 








fti(0) 





h (D 


fto(O) 




hi(l) 


hi(0) 




fto(l) 







feid) 


h (L h - 1) 




fto(0) 


: ^(L h - 1) 







ho(L h -i) \ 




ho(l) 





feiaji - 1) 





h (L h - 1) 






hi(0) 
hid) 

hi(L ft - 1) 



then (2.21) can be written as below 



(2.23) 



w(n) = (H^g}w(n) = [H' i H'J 



go 
gi 



•w(n). 



(2.24) 



w(n) = < 



fto(O) 
h (l) 



fto(O) 





ha(0) 
hid) 



hxCO) 


fto(L ft -l) 






ho(l) 

ho(i ft -l) 




h (0) 

i fto(l) 

... ft (L h -l) 


hiO*-!) 






hid) 

hxCLft-1) 







hx(0) 

hxCLft-1) 



r 00 (o) pi 


5od) 


5o(i 5 -l) 


0i(O) 


0i (1) 


^i(L fl -l)l 



• w(n). 



(2.25) 
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with 



w 



(n) = [w(n) w(n - 1) ... w(n -L h -L g + l)f . (2.26) 



where H T and H T 1 are the convolution matrices of the two subchannels. Note that 
there is no need to know the channel order L h precisely, only the equalizer length L g 
has to be chosen so that the following hypotheses hold in order to have sufficient 
channel diversity 

H t ) L h + L g — 1 < PL g (H T has strictly more columns than rows). 

H 2 ) H T has full row rank L h + L g — 1 which means the matrix H T is invertible. 

The length of equalizer L g must be chosen large enough (i.e. L g > L h — 1) in 
order to have sufficient equations to solve the equalizer. In this case, a perfect 
equalization can be achieved. It has been shown as in [8], [10], [11], [18], and [24] 
that the block Toeplitz structure of H T implies that it will be full row rank if and only 
if there are no common subchannel roots and L g > L h — 1. It is also found that there 
is a unique solution of the equalizer equations if L g = L h — 1 . 
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CHAPTER III 
BLIND ADAPTIVE CHANNEL EQUALIZATION METHODS 

3.1 Adaptive Filter Theory 

(The beginning of this section is inspired from [13].) 

An adaptive filter is defined as a self-designing system that relies for its 
operation on a recursive algorithm, which makes it possible for the filter to perform 
satisfactorily in an environment where knowledge of the relevant statistics is not 
available. Linear and non linear are main groups of adaptive filters. Linear adaptive 
filters compute an estimate of a desired response by using a linear combination of the 
available set of observables applied to the input of the filter. Otherwise, the adaptive 
filter is said to be non linear. In this chapter, adaptive filters theories are described for 
different kinds of algorithms. 

Adaptive filters can be integrated in systems with different functionalities, 
being predictive deconvolution, system identification, channel equalization, spectral 
analysis, signal detection, beamforming, interference and noise cancellation examples 
of such. Also, adaptive filters can adjust to unknown environment, and even track 
signal or system characteristics varying over time. 



22 



3.2 Adaptive Equalization 

In digital communication systems, one application of adaptive filters is 
adaptive channel equalization that removes ISI generated from unknown channel. ISI 
can be compensated or removed from received signals by using many approaches 
such as optimum receiver which uses maximum likelihood estimation, suboptimum 
receiver which uses a linear equalizer or a decision-feedback equalizer [20]. The three 
equalization methods assume impulse response of the channel characteristics or 
frequency response of the channel characteristics were known a priori to the receiver. 
However, in most practical digital communication systems that use equalizers the 
channel characteristics are not known a priori. Therefore, training sequences that 
provide different realizations of a desired reference signal can be used to estimate the 
channel and find the necessary equalizer. The reliance of an adaptive channel 
equalizer on a training sequence requires that the transmitter cooperates by resending 
the training sequence, lowering the effective data rate of the communication link [7] . 

In high speed digital communication systems, the transmission of a training 
sequence is either impractical or very costly in terms of data throughput. Another 
drawback of using training sequences for channel equalization is that time and 
bandwidth are consumed for the equalization process, causing most adaptive filters 
whose desired reference signals are training sequences cannot be used. For these 
reasons, blind adaptive channel equalization algorithms have been invented and 
developed, and they have been the best approaches that satisfy our system models. 
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The purpose of this thesis is applying blind adaptive channel equalization with the 
variety of algorithms. 

3.3 Blind Adaptive Equalization 

In some applications, such as multipoint communication networks, the input 
signal to the channel is unknown to the receiver. Only statistical properties of the 
input structure that blind equalization based on to recover the unknown input 
sequence are known. Therefore, it is desirable for the receiver to synchronize the 
received signal and to adjust the equalizer (self-adaptation) without having a known 
training sequence available. Equalization techniques based on the initial adjustments 
of the coefficients without benefit of the training sequence are said to be blind. Now, 
we can define a blind adaptive filter is an equalizer that uses one of these techniques 
to perform recursive adjustments of its parameters without the aid of training 
sequences. 

3.4 Blind SIMO Channel Equalization Methods 
Using Second Order Statistics (SOS) 

In the blind methods of equalization, some structure properties of the input 

signal are used such as whiteness in conjunction with the receiver outputs in order to 

estimate the equalizer. Also, most commonly used adaptive algorithms for blind 

channel equalization do not require extra bandwidth for training. Therefore, blind 

equalizers have gotten great research and practical interest, and many effective 

algorithms have been proposed for last few years [19] and [23]. In recent papers, ( see 

for example [3], [18], [26] and [27]), it has been shown that the second order statistics 
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(SOS) of the channel output y(ri) is sufficient to estimate most communication 
channels when the received signals are fractionally sampled (oversampling), or 
multiple antennas are used. 

The methods introduced in [9], [10], [11] and [22] have been used SOS of the 
channel output to directly estimate the equalizer without finding the channel impulse 
response as a first step. These types of methods are called blind methods. In the 
following section, the cyclic statistical property of the fractionally sampled 
observation [y(n)] is described because some calculation algorithms that will be 
discussed later in this thesis use this property. 

3.5 Cyclic Statistics 

In this section, we describe the co variance of y N (n) which will be useful for 
finding the equalizers directly from the noise-free channel output. The following 
vectors are defined to represent N vector observations of the SIMO fractionally 
spaced model (see Figure 2.5). 

yjvOO = [y 00 y (n - 1) ... y Q (n - N + 1) y^n) y x {n - 1) 

... yi (n-N + l)Y, 
*w( n ) = [ x o( n ) x o( n ~~ 1) ■■■ x o( n — N + i) x-tCn) x-^n — 1) 

... Xi(n - N + 1)] T , 
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^jv( n ) = [^o( n ) v o(. n — 1) ■■■ ^o( n — N + 1) v^n) v x (n — 1) 

... v^n-N + 1)] 7 . (3.1) 

The correlation of the scalar output y(n) in equation (2.6) is defined 

c 2y (n; m) = E{y(n)y*(n + m)}. (3.2) 



oo oo 



c 2y {n; rri) = y } c 2w {l 2 - lj)h(n - 2l 1 )h*{n - 2Z 2 ) 

l 1 = — co l 2 = — co 

+c 2v (ni). (3.3) 

where m = 0,1, ...,L g — 1, c 2w (m) = E{w(n)w*(n + m)} and c 2v (v) = 
£"{i;(n)v*(n + m)} are the correlations of the stationary input and noise, respectively. 
It can be proven from (3.3) that the correlation is periodically time-varying in n with 
period 2, see for proving [25] . 

c 2y (n; m) = c 2y (n + 21; m) V I integer. (3.4) 

Due to the periodicity of the correlation matrix, we can use y N (n) in similar 
manner instead of using y(n) to find the correlation matrix C 2y as below 

C% = E{y N {n)y N {n) H }. (3.5) 
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3.6 Direct Blind MMSE Equalizers 

In this section, we address noise suppression and consider FIR Weiner filters 
to determine the minimum mean-square sense estimate of equalizer output w(n) 
using only y(n). The objective is to compute the equalizer g(n) so that the following 
cost function/ is minimized. 

/ = E{\w(n) - w(n - d)\ 2 }. (3.6) 

where d is the integer delay between input w(ri) and noise-free channel outputs 
Xj(n), i = 0,1- The output of equalizer w(n) can be written as below 



Hn) = tC g Q)yo(n - + ti' 1 giQ)yi(n - 1), 



w(n) = Y.Ut.C 9iH)yi{n - I). (3.7) 

We substitute (3.7) into (3.6) and take the first complex partial derivative with 
respect to the unknown equalizer coefficients and set it equal zero to minimize / 



a 



E^Ut.C gmytin - I) -win - d)| 2 } = 0. (3.8) 



dg* k {m)' 

where k = 0, 1 and m = 0,1, ...,L g — 1. After simplification (see [25]), it yields to 
orthogonality condition 



1 L 3 1 



^ ^ g t (i)E{ yi {n - 0y k *(n - m)} 



i=0 1=0 

-E{w(n-d)y k *(n-m)} = 0. (3.9) 
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Since the input symbols are assumed to be i.i.d., and uncorrelated with the 
noise, the second term of (3.9) is exist if and only if I = n — d. Therefore, the second 
term can be written as follow 



h ~ 1 
E{w(n - d)y k *(n - m)} = E \ w(n - d) > w*(l)h* k (n -m-l)>, 

1=0 



- E{w(n - d)w*(n - d) 

x h* k (n — m — (n — d))}, 

= E{w(n - d)w*(n - d)}h* k (d-m), 

= alh* k (d-m). (3.10) 

substituting (3.10) into (3.9), yields 

l L a~ 1 
2^ 2^ g t (l)E{yi(n - l)y k *(n - m)} - <r£h* k (d - m) = 0, 

i=0 1=0 

lUlilo 1 Si(0E{y t (n - l)y k *(n - m)} = <r*h' k (d - m). (3.11) 

the above equation can be written in a vector form as 

E{y* Lg (ri)yl g (ri)}g d = <tfH*(:,d + 1). (3.12) 

C 2 *?g d = <tfH*(:,d + l). (3.13) 

go=^[C 2 *y] + H*(:,l). (3.14) 
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ga = ti[C?y?C£ d [C 2 ^]tH*(:,l). (3.15) 

The above equations (3.14) and (3.15) give to us the zero-delay MMSE 
equalizer by substituting d — and the arbitrary nonzero delay (d =£ 0) MMSE 
equalizer. 

3.7 Methods for Blind Equalizer Calculation 

This section introduces and describes five methods for blind equalizer 
calculation with their derivations. These blind algorithms are procedures for self 
adjusting the parameters of an adaptive filter to minimize the cost function. We 
describe the general form of adaptive FIR and IIR filtering concept. 

3.7.1 Recursive Blind Adaptive Using Recursive Least 
Square Algorithm (RLS) 

This algorithm addresses the problem of blind adaptive equalization of finite 
impulse response channels (FIR) with the exploiting the diversity induced by sensors 
arrays or fractionally sampled. It analyses and solves the case when we have parallel 
subchannels in the presence or absence of noise. It is shown later that the resulting a 
priori error converges towards a scalar version of the input symbol sequence and the 
equalizer parameters are estimated by using Recursive Least Square (RLS) algorithm. 
This algorithm is very useful for tracking time varying channel. The RLS algorithm 
uses information from all past input samples (and not only from the current tap-input 
samples) to estimate the static correlation matrix of the input vector C^y ■ To start 
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driving this algorithm, it is assumed that the additive noise is to be uncorrelated with 
input samples and a forgetting factor A (0 < A < 1) is included to reduce the effect of 
past observations on the statistics correlation matrix estimate C^y ■ 

c\ y {N) = E{yl g {n)y T Lg {n)l 



1=0 



At / = N, the correlation matrix is simplified as 



N-l 



C 2y (N) = £ A N - l yZ g (l)y T Lg (l) +yt g (N)y T Lg (N) 



1=0 



= AC* y (N - 1) + y* Lg (N)yl g (N). (3.16) 

The result of (3.16) is a recursive form, and it can be applied to (3.14) to find 
MMSE equalizer g with zero delay. However, using the matrix inversion lemma [13, 
p. 480], a direct method to find g is developed without requiring a matrix inverse or 



-l 
and 



pseudo in version of C^y ■ This can be done by defining P(iV) = C^CiV) 
applying the matrix inversion lemma with (3.16) 

P{N) = A^PiN - 1) 

A~ 2 P{N - i)y* L (N - DyllN - 1)P(JV - 1) 

(3 17^) 

1 + A-^yI g (N - 1)P(JV - l)yt g (N - 1) ' K ' } 
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The above equation can be applied to (3.14) to estimate the MMSE zero delay 

A 

equalizer g (N) 



g (JV) = POV)iA_ ff 2H r (l,:). (3-18) 



The equalizer estimate form in (3.18) does not require a matrix inverse and, 
hence, is computationally feasible for adaptive implementation. If A = 1 is chosen, 
(3.18) provides a method to recursively compute the time invariant equalizer taps, 
thereby reducing the memory requirements for long data records [11]. It is very 

A 

important to mention that the convergence of C 2y (iV) is affected by the forgetting 
factor A and, thereby, it affects the accuracy of the estimator. The forgetting factor A 

A 

is usually chosen between [0.98, 1]. The delay MMSE equalizer gd(iV) can be driven 
in a similar way using the same above steps in this method. 

3.7.2 Recursive Blind Adaptive Using Cyclic Least 
Mean Square Algorithm (cyclic LMS) 

This algorithm also considers the problem of blind adaptive equalization of 

finite impulse response channels (FIR) with the exploiting the diversity induced by 

sensors arrays or fractionally sampled. It analyses and solves the case when we have 

parallel subchannels in the presence or absence of noise. The cyclic LMS algorithm is 

very popular, and it uses the Stochastic Gradient Descent approach for updating the 

equalizer coefficients at each symbol. It provides an instantaneous approximation 

estimate to the gradient vector of the cost function. The goal of the cyclic LMS 

method is to estimate the equalizer coefficients to achieve the least mean squared 
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error. It is shown later that the resulting error converges towards a scalar version of 
the input symbol sequence. The update of the coefficients is performed in the 
following gradient descent equation 

goOV) = g (JV - 1) - ^V/oOV). (3-19) 

where u is the step size parameter which controls the moving distance along the error 

A 

surface, and VJ (N) is the instantaneous approximation to the gradient of the cost 
function J = E{|w(n) — w(n)| 2 }. 



A 



V/ = 



dJo dj dj dj dj dj 



La g *(o) a g j(o) agS(i) agj(i) agSa^-i) d % \(_i g -i)\ 



(3.20) 



The equation (3.20) can be simplified in a same manner as in the derivation of the 
MMSE equalizer (see section 3.6), yields 

V/o =E[yt a yf„}go-<tfH"(l,:). (3.21) 

as a result, the instantaneous approximation at time JV is obtained by 



V/o(JV) = y* Lg Wyl g (N)g (N - 1) - <^ H"(l, 0- (3-22) 



substituting (3.22) into (3.19), yields to the following equation that the equalizer 
coefficients can be estimated 

goOV) = go(JV - 1) - ln\yt g (N)yl g (N)g (N - 1) - <r*H H (l, :)](3.23) 
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The speed of convergence and steady-state performance rely on choosing of 
the step size \i. The convergence analysis issues are not addressed in this thesis, and 
they have been studied in many papers and books such as [13]. 

The cyclic LMS algorithm, as we can see through derivation, does not rely on 
an explicit matrix inverse. Therefore, it is not sensitive to nearly common subchannel 
roots. In addition, it has extremely low computational complexity, but it has slow 
convergence [13, pp. 334- 335]. 

3.7.3 Recursive Blind Adaptive Using Least Mean 
Square Algorithm (LMS) 

The LMS algorithm is widely used in various applications of adaptive filtering 
due to its computational simplicity. This algorithm is similar to previous section, the 
step size parameter fi, which affects the convergence speed of the LMS, is used to 
control the moving distance along the error surface. They use information only from 
the current tap-input symbols. Also, they both depend on the statistics of the input 
w(n) and the output y(n) signals since their updating equations consist of first 
complex partial derivative to the cost function that has expectation to those two 
signals. However, this algorithm uses a Steepest-Descent-Based algorithm for 
updating the equalizer coefficients at each symbol [9] . 

In section 3.6, the optimal solution for parameters of the adaptive filter is 
driven. The objective of this optimal (Wiener) solution is to compute the equalizer 
g(n) by minimizing the cost function/ relying only on using the observation y(n). In 
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other word, this solution leads to the minimum mean-square error in estimating the 
reference signal w(ri). The optimal (Wiener) solution is given by 

go = [C^E { W {n)y\n)l 

go = tfw[C2y] +H *C,l). (3-24) 

where C^y = E{y^(ri)yJi(n)} and let p = E{w(n)y^(n)} which both are unknown. 
Now, the LMS method updates the equalizer coefficients using a steepest-descent- 
based algorithm which in turn can be used to search the Wiener solution of equation 
(3.24) as follows 



A 

g(n + l) = g(n)-nJ(n), 



g(n + 1) = g(n) + 2n[p(ri) - C$ (n)g(n)]. (3.25) 



A 



where J(ri) is an instantaneous estimate to the gradient vector of the cost function 
with respect to the filter coefficients. Employing instantaneous approximation 



A 



estimates C 2 *y ( n ) and p(n) for C^y and p is one possible solution to estimate the 
gradient vector as follow 



C${n)=y* N (n)y T N {n). (3.26) 



A 



p(n) = w(n)y* N (ri). (3.27) 
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The gradient estimate is given by 



/(n) = -2[p(n)-C 2 *J r (n)g(n)] 



= -2w(n)y* N (n) + 2y* N (n)yl{n)g{n) 

= -2y* N (n)[w(n) -y^(n)g(n)] 

= -2y^(n)err(n). (3.28) 

This can be applied back to (3.25), the resulting is the updating equalizer parameters 
estimate of the least mean square LMS algorithm. In summary, the updating 
equations for the LMS algorithm are described by 

err(n) = w(n) — yJi(ri)g(n). (3.29) 

g(n + 1) = g(n) + 2ny* N (n)err(n). (3.30) 

The LMS coefficients update which is illustrated in the equation (3.30) is a 
form of time-averaging that smooth the errors in the instantaneous gradient 
calculation to obtain a more reasonable estimate of the true gradient [9]. 
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3.7.4 Unbiased Blind Adaptive Using Gradient Projection Technique 

This algorithm considers the problem of blind adaptive equalization of finite 
impulse response channels (FIR) with the exploiting the diversity induced by sensors 
arrays or fractionally sampled. It analyses and solves the case when we have parallel 
subchannels in a noisy environment. An adaptive equalizer is developed depending 
on the FIR quadratically constrained filter which uses the Gradient Projection (GP) 
technique as a solution for its implementation at the receiver. It has been shown that 
the (GP) technique as described for example in the work by [6] is the most direct 
solution to the adaptive implementation of quadratically constrained filters. Also, the 
Gradient Projection algorithm can be used in adaptive beamforming problems at low 
cost sense [10]. It is shown later that the resulting prediction error converges towards 
a scalar version of the input symbol sequence and the adaptive filter parameters are 
estimated by using (GP) algorithm. In addition, the additive noise is assumed to be 
uncorrelated with input samples. 

In order to start driving this method, recall form chapter two that there is no 
need to know the channel order L h precisely, only the equalizer length L g has to be 
chosen such that the following hypotheses hold in order to have sufficient channel 
diversity 

H ± ) L h + L g — 1 < PL g (H T has strictly more columns than rows). 

H 2 ) H T has full row rank L h + L q — 1 which means the matrix H T is invertible. 
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The notion of linear prediction method (LP) is whitening the observation to 
find a channel inverse filter, using the following a priori constraint: 

H 3 ) The input sequence w(ri) is white, E{w(n)w*(n — k)} = 8 k . 

It is a good point to mention that there is limitation in the SISO case because 
the whiteness constraint is too weak to allow phase and amplitude equalization of a 
mixed-phase channel. However, in the SIMO model, this limitation vanishes due to 
hypotheses H t and H 2 , so the exact FIR inverse of the nonideal channel exists. It has 
been shown in [10] and [24] that the resulting noise-free prediction error signal 
related to the input symbol sequence w(ri). 

Assume that the finite-length equalizer vector g be a 2L g x 1 complex- valued 
vector of prediction coefficients, and err k (ri) be the prediction error defined for k = 
0,1- The prediction error achieves its optimality in the noise free case if and only 
if err k (n) = h k (0)w(ri), and thus the prediction error variance /(g) = E\err k (ri)\ 2 
is minimized. The performance versus SNR of this theoretical equalization method 
would critically depend on the realization of the particular coefficient h k (0), so it is 
more appropriate to exploit the predictor as a tool to identify the channel [10]. A 
major drawback of the prediction error filtering is sensitivity to channel additive 
noise. For this reason, a modified prediction scheme will be introduced to allow the 
adaptive computation of an unbiased predictor in a noisy environment. Adaptive filter 
parameters are estimated by using the Gradient Projection technique. The algorithm 
which is first introduced by [10] goes as follows: 
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err k (n) = y k {n) - g(n - l) H y N (n - 1). (3.31) 

UfcOO =y N (n- l)err k (ny. (3.32) 

g(n) = g(n - 1) + au k (n). (3.33) 



*» - KB) VlK-)l- ( » 4) 



where a is a small step size and 7 is a constraint value. The performance behavior of 
the mean square error is influenced by the choice of y in this particular algorithm. It 
is not difficult to see that (3.33) is a standard LMS update of equalizer parameters. 
This algorithm has a very low computational cost, and it shows desirable robustness 
properties [10]. 

3.7.5 Recursive Blind Adaptive Using Recursive Extended 
Least Square Algorithm (RELS) 

This algorithm considers the problem of blind adaptive equalization of infinite 

impulse response (HR) channels without requiring the channel diversity condition, 

and it analyses and solves the case when we have parallel subchannels having 

common zeros. This common factor is assumed as a minimum phase filter, while 

overall subchannels can be a non minimum phase systems. An equalizer is developed 

based on the optimal IIR filter as a predictor of the received signal. The main criterion 

of this method is a one- step ahead prediction of one of the subchannel outputs [1], 

[2], [10] and [22]. It is shown that the resulting prediction error converges towards a 

scalar version of the input symbol sequence and the adaptive predictor parameters are 
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estimated by using Recursive Extended Least Squares (RELS) Algorithm. In addition, 
it is assumed that the input samples are mutually uncorrelated, and the parameter 
estimates are updated when each single signal is received [22] . 

Prediction -Based Equalizer 

Consider the case of single-input two-output system model which means the 
receiver performs two measurements for each transmitted symbol. An equivalent 
representation to our channel model is shown in Figure 3.1. 



Xl(l) = mt*> w(0 ' 



*2(0 = A f -i-v w(0- (3-35) 



where q 1 is a unit delay operator, w(i) is the transmitted symbol and 
i4(q _1 ) =l+a 1 q~ 1 + — \- a L q~ L , 

Biq-^bo+^q- 1 + .» + b L q~ L , (3.36) 

C(q _1 ) =Co+Ciqf _1 + ••• + c L q" L , 

F(q- 1 )=i+/ 1 <r 1 + -+/nF<r nF . 

where L is the polynomial channel order. Assuming that polynomials i4(q _1 ), 
5(q _1 ), and C(q _1 ) are of the same order, and assuming that / A( n -i\ is a stable 
operator. Also, S(q _1 ) and C(q _1 ) are coprime polynomials as well as F(q _1 ) is a 
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minimum phase polynomial. In general, all quantities in (3.35) and (3.36) can be 
complex numbers. 



w(i) 





Bdq-^Fiq- 1 ) 


*i(0 






A(q-i) 




















Ciq-^Fiq- 1 ) 
A(q-i) 


x 2 (i) 









Figure 3.1 SIMO channel model (IIR). 
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Figure 3.2 Predictor based equalizer. 
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The prediction-based equalizer is shown in Figure 3.2, where y(i + 1) is one 
step-ahead prediction of x 1 (i + 1). Let i?(q _1 ), 5(q _1 ) and D(q _1 ) be the filter 
operators and the polynomials in q _1 , so y{i + 1) is 



y(l + 1} = D(^j * l(0 + D(^ X2( °- (337) 



which is optimal in the mean-square sense. The polynomials R(q x ), 5(q x ) and 
D(q _1 ) are computed by minimizing the following cost function/ : 

J = E{\x 1 {i + l)-y{i + l)\ 2 ). (3.38) 

Note that we can use x 2 (i) instead of x 1 (i) as a reference signal and find the 

predictor by minimizing £"(|x 2 (i + 1) — y(i + 1)| 2 ). From (3.35) and (3.37), 
x t {i + 1) can be written in the following form 

x x (i + 1) - y(i + 1) = -x(i) + b w(i + 1). (3.39) 



where 



x(i) = H(qf~ 1 )w(0- (3-40) 



and 



Aiq- 1 ) 



(3.41) 
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since F(q _1 ) is not a factor of i4(q _1 ), from (3.41) it follows that #(q _1 )=0 if and 
only if D(q _1 ) = F(q _1 ) whereas i?(q _1 ) and 5(q _1 ) are solutions of the following 
polynomial equation: 

[BCq-^RCq-^ + CCq-^SCq- 1 )] 

= qtfdq-^FCq- 1 ) - ^(q" 1 )^] (3-42) 

notice also that H(q _1 ) = gives x(i) = and from (3.39), we get the prediction 
error 

Xi(i + 1) - y(i + 1) = 6 w(i + 1). (3.43) 

Let the predictor operators i?(q _1 ) and 5(q _1 ) are defined as 

i?(q _1 ) = r + r^" 1 + ••• + r Ni q~ N \ 

Siq- 1 ) =s + s^- 1 + ... + s^qf" 2 . (3.44) 

Then (3.42) has a solution with respect to i?(q _1 ) and S(q _1 ) if A^ > JV , 
iV 2 > iV , where N Q = max(L — l,nF — 1). A unique solution exists if 5(q _1 ) and 
C(q _1 ) are coprime, L > nF and at least one of the following holds: N-^ = 

N or N 2 = N . 

In practice, we cannot use (3.42) to calculate i?(q _1 ), 5(q _1 ) andD(q _1 ) 
since polynomials 5(q _1 ), C(q _1 ) and F(q _1 ) are unknown. Therefore, we cannot 
use a nonadaptive predictor (3.37) to calculate prediction error (3.43). A recursive 
algorithm is proposed for directly estimating the unknown parameters in (3.37). 
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Taking advantage to our assumption that D(q x ) = F(q 1 ), the nonadaptive 
predictor (3.37) can be written in the form 

Fiq-^yii + 1) = R(q- 1 )x 1 (0 + S(q~ 1 )x 2 (i)- (3-45) 

Let 

0o(O T = [*i(0< -,*i(i - N 1 ),x 2 {i), 

...,x 2 {i-N 2 ),-y{i) -y{i-N 3 )l N 3 >nF-l, (3.46) 

and 

"o = L r o< ■■■' r w 1 ' 5 0' ■■■' s N 2 'h' —ifnFi 0' ■■■» 0J. (3.47) 

Then (3.45) becomes 

y(i + l)= OoCO- (3.48) 

The number of inserted zeros in (3.47) relies on the value of JV 3 and the 
unknown degree nF. Equation (3.48) is called the optimal nonadaptive prediction 
because 9 is unknown. Therefore, the following adaptive predictor can be used 

y(i + 1) = §(i + l) H 0(i)- (3.49) 

where 9(i) is an estimate of unknown 9 , and 

0(i) T = [Xi(0- ...,x 1 (i-iV 1 ),x 2 (0, 

...,x 2 (i-N 2 ),-y(i) -Ki-N 3 )]. (3.50) 
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Equation (3.49) is called a posteriori prediction, whereas y{i + 1) = 9(V) H 0(i) 
represents a priori prediction. A priori prediction error x t {i + 1) — y{i + 1) is used 
to run the following extended least squares algorithm (first presented by [22]): 

0(i + 1) = 0(0 + p(i)0(i)e(i + 1)*, (3.51) 

e(i + 1) = Xi(i + 1) - 0(O H 0(O, (3.52) 



p(l) = p(l - 1} - i + 0(O» P G- D0(O ■ (3 " 53) 



p(0) = p I, Po > 0, 

where I is the identity matrix, and p is an arbitrary finite positive scalar. Initial 
values of 0(0) in (3.51) and y(k), k < in (3.50) are chosen arbitrary. Predictor 
coefficients are directly computed using the RELS algorithm (3.49)-(3.53). There is 
no need to know the channel order L and nF precisely, as long as one overfit, i.e., in 
(3.50), JVi > N , N 2 > N , N = max(L - 1, nF - 1), and JV 3 > nF [22]. 
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CHAPTER IV 



PERFORMANCE OF ADAPTIVE BLIND EQUALIZATION 
ALGORITHMS AND SIMULATIONS 



This section presents the results of simulations using MATLAB to examine 
the performance behaviors of various adaptive /recursive algorithms described in 
chapter 3. The performance of adaptive blind algorithms is assessed through this 
thesis by calculating the mean-square-error (MSE), and showing constellation plots of 
each single method. The major mean of comparison is the error cancellation 
capability using algorithms that rely on the parameters such as step size fi, forgetting 
factor X, number of iterations and their MSE performance. The principle advantage of 
these algorithms is to remove ISI, which is generated during transmission, by 
equalizing the transmission channels, or we can say in a more fashionable way it is to 
reduce the error probability in the decision at the receiver. The calculation methods 
that will be simulated in this chapter are cyclic LMS, LMS, RLS, Linear Prediction 
(using GP technique) and RELS. Before the simulation results are presented, it is 
useful to look at the brief explanation to the MSE and symbol constellations as 
performance measures in the following section. 
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4.1 Performance Evaluations 

4.1.1 Mean Square Error (MSE) 

In general, the minimum mean square error (MSE) is a metric indicating how 
well a system can adapt to a given solution. A small minimum MSE is an indication 
that the adaptive system has accurately modeled, predicted, adapted, and/or 
converged to a solution for the system. In other words, the algorithms will achieve 
better performance. A very large MSE usually indicates that the adaptive filter cannot 
accurately model the given system, or the initial state of the adaptive filter is an 
inadequate starting point to cause the adaptive filter to converge. The ideal MSE is 
when it reaches zero. However, most real systems cannot achieve the ideality. The 
MSE is defined in a similar form as a cost function in (3.6) 

MSE = E{\w{n-d)-w(n)\ 2 }. (4.1) 

This equation can be simplified and used in practice by a consistent sample 
estimate based on JV observations 

MSE{N) = ^lSlw(n - d) - w(n)| 2 . (4.2) 

where d is the desired delay. Most of the time, the graphs illustrate the MSE versus 
time index samples. The MSE achieves its asymptotic steady state level after a 
sufficient number of symbols. 
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4.1.2 Symbols Constellation Plots 

It is also called scatter or eye diagram. This type of performance measures is 
usually useful to see if the algorithm used is working or not. It is a constellation 
diagram of the observed signals and equalizers outputs. 
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4.2 Simulations 

In this section, the simulation results are presented for all adaptive /recursive 
methods that were discussed in the previous chapter. The calculation methods that 
include cyclic LMS, LMS, RLS, LP and RELS are simulated using MATLAB. Some 
experiments assumed that the input signal is independent and identically distributed 
(i.i.d.) sequence with the zero mean and unit variance; however, other experiments do 
not exploit this identity. Nine experiments will be performed in the following nine 
sections such that each section includes graphs and results related to each single 
calculation method. 

4.2.1 Experiment la 

In this experiment, an i.i.d. symbol sequence that is generated from a 16 QAM 
constellation is used. The symbol levels along both axes are -1.5, -0.5, 0.5 and 1.5. 
The continuous time channel that will be the first class of channel model used in the 
simulation of this thesis spans four symbols and describes for t G [0, 4T S ) 

h c {t) = e-^(° 15 )r c (t - 0.25r s ,/?) + 0.8e^' 27T ^ a6 )r c (t - T s ,p). (4.3) 
where r c (t,/3) is the raised cosine given in [20, pp.546] as 

wt COS (— -) 

r c {t, /?) = sinc(-) — ^Hr- (4.4) 

cv >rj y T s J l-4:l3 2 t 2 /Tl v ' 

with roll-off factor /?, while T s is the symbol duration. As in [11], /? = 0.35 is 
chosen. 
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The above h c {t) represents the composite causal approximation of a two-ray 
multipath mobile radio environment. This channel is reported by [1] and [21]. The 
discrete-time equivalent channel h(ri) is obtained by oversampling (P = 2) h c {t) at a 
rate of T s /2 (fractionally sampled) or h(n) = h c (nT s /2) for n = 0, 1, ...,7. Figure 
4.1 shows the magnitude of the impulse response h(n) while Figure 4.2 shows the 
zeros of h(ri), and the zeros of the subchannels h (n)and /i 1 (n). The subchannels 
parameters are 

h = [0.52 -jO.72 - 0.48 +;0.24 - 0.05 +;0.07 0.01 -;0.02] r . 

(4.5) 

h x = [0.12 -yo.43 - 0.48 +;0.41 0.13 -7O.il - 0.04 +;0.03] T . 



(4.6) 
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Figure 4.1 Two-ray multipath channel, the magnitude of impulse response h(n). 
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Figure 4.2 The zeros of h(ri), and the subchannels zeros h (ri)and /i 1 (n). 

Both Figures 4.1 and 4.2 are first presented in [11]. Using this type of channel 
to implement the first equalizer calculation method that is Recursive Least Square 
(RLS) with zero delay (d = 0). It is assumed that the subchannel length L h — 4, and 
subequalizer length L g = 4. The received constellation is presented in Figure 4.3, 
while Figure 4.4 shows the equalized eye diagram. Also, Figure 4.5 depicts the mean- 
square symbol error for RLS algorithm across 3000 symbols. All of these figures in 
the presence of the receiver noise v(ri) that is assumed to be white with zero mean 
and 0.14 variance. The performance of RLS algorithm depends on the forgetting 
factor X which in this experiment is chosen to be 0.99. For simulation, use equations 
(3.17) and (3.18). 
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eye diagram of y1 channel output 




Figure 4.3 Constellation plot of the received sequence for 16 QAM 

modulation. 
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Figure 4.4 Equalized symbol-eye diagram, using RLS algorithm with the raised- 
cosine channel. 
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Figure 4.5 Mean-square symbol error of RLS algorithm. 

We can see from Figure 4.5 that this algorithm does a very good job of 
equalization after as few as 1000 symbols. The MSE converges to 0.06665 at index 
sample n = 3000. 

4.2.2 Experiment lb 

In this section, an uncorrelated symbol sequence generated from a 16- QAM 
constellation is used. The channel model (second class), which is used in this 
experiment, is reported by [22], and it has different coefficients than the previous 
experiment. Figure 4.6 shows the magnitude and zeros of the impulse response h(ri). 
The subchannels parameters are described as follow: 
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h = [0.5498 -y0.3174 0.4342 -;0.4596 0.0131 -yO.0919 



-0.0275+70.1592]' 



(4.7) 



h x = [0.0414 - ;0.0717 0.7168 - ;0.3407 0.0745 - ;0.3636 



0.0685 -;0.1861] 5 



(4.8) 
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Figure 4.6 The magnitude and zeros of impulse response h(ri), second class of 

channel model. 



The RLS algorithm is implemented again using the second class of channel 
model with zero delay (d = 0). It is assumed that the subchannel length L h = 4, and 
subequalizer length L g = 4. The received symbols are presented in Figure 4.7, while 
Figure 4.8 shows the equalized constellation. Also, the mean-square symbol error is 
estimated for this algorithm across 3000 symbols as shown in Figure 4.9. All of these 
figures in the presence of the additive noise v(n) that is assumed to be white with 
zero mean and 0.14 variance. The forgetting factor X is also 0.99. 
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Figure 4.7 Eye diagram of the received signal for 16 QAM modulation. 
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Figure 4.8 Equalizer output symbols, using RLS algorithm with the second class of 

channel model. 
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Figure 4.9 Mean- square symbol error of RLS algorithm. 

From Figure 4.9, we can see that this algorithm works well when the second 
class of channel model is used. The MSE converges to 0.0709 at index sample 
n = 3000. 

4.2.3 Experiment 2 

In this experiment, the cyclic LMS algorithm, which is first presented in [11], 
is simulated with zero delay (d = 0) using the first class of channel model (the 
raised-cosine two-ray multipath mobile radio) with L g = L h = 4. Figure 4.10 shows 
the equalized eye diagram, whereas Figure 4.11 depicts the mean-square symbol error 
for cyclic LMS algorithm across 3000 symbols. The additive noise is the same as last 
experiment. Trading off speed convergence with steady-state error, the performance 
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of cyclic LMS relies on the step size \i which throughout this thesis is chosen to be 
0.0025. For simulation, use equation (3.23). 
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Figure 4.10 Eye diagram of equalizer output, using cyclic LMS with the raised-cosine 

channel. 
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Figure 4.11 Mean-square symbol error of cyclic LMS. 
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From Figure 4.11, the performance of cyclic LMS is better than RLS 
algorithm since cyclic LMS has lower MSE and faster convergence. The MSE 
converges to 0.02526 at index sample n = 3000. However, the cyclic LMS doesn't 
work with the second class of channel model. To make it work, the conjugate 
calculation y[ (iV) in (3.23) has to be removed. Otherwise, it won't work. 

4.2.4 Experiment 3a 

In this experiment, The LMS algorithm is performed using the raised-cosine 
channel with zero delay (d = 0). We still suppose the same assumption as in the 
experiment 2. Figure 4.12 shows the equalized constellation, while the mean-square 
symbol error is estimated for this algorithm across 3000 symbols as shown in Figure 
4.13. For simulation, use equations (3.29) and (3.30). 
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Figure 4. 12 Equalized eye diagram, using LMS algorithm with the raised-cosine channel. 
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Figure 4.13 Mean-square symbol error of LMS algorithm. 

As we see from Figure 4.13, the LMS method works well. The MSE 
converges to 0.05596 at index sample n = 3000. 

4.2.5 Experiment 3b 

In this section, the second class of channel model is used to implement the 
LMS algorithm. Our assumptions are still the same as in previous section. Figure 4.14 
shows the equalized eye diagram, whereas Figure 4.15 depicts the mean-square 
symbol error for LMS algorithm across 3000 symbols. 
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Figure 4.14 Eye diagram of equalizer output, using LMS algorithm with the second 

class of channel model. 
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Figure 4.15 Mean-square symbol error of LMS algorithm. 
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From Figure 4.15, the performance of this experiment is better than 
experiment 3a since the MSE of LMS using the second class of channel model has 
lower value at 3000 (index sample). For more clarity, the MSE converges to 0.04817 
at index sample n = 3000. 

The comparison between RLS, cyclic LMS and LMS is accomplished through 
Figure 4.16 that shows the MMSE for the cyclic LMS, LMS and RLS. The channel 
used for this comparison is the raised-cosine two ray multipath mobile radio (P = 2). 
The plot demonstrates that the cyclic LMS has better performance since it has the 
lowest MSE. The MSE of LMS is lower than the MSE line of the RLS at (n = 3000). 
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Figure 4.16 Mean-square symbol error, Comparison using raised-cosine channel. 
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The second comparison is also between cyclic LMS, LMS and RLS, but here 
the second class of channel model is used. We can see from figure (4.17) that the RLS 
has the best performance due to its MSE, which is lower than the MSE of both cyclic 
LMS and LMS. As I mentioned earlier that the cyclic LMS doesn't work with this 
channel class and it will not converge. 
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Figure 4.17 Mean-square symbol error, Comparison using second class of channel 

model. 
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4.2.6 Experiment 4a 

In this experiment, the unbiased blind adaptive using Gradient Projection (GP) 
technique is implemented. It is first presented in [10]. The channel used is the raised- 
cosine two-ray multipath mobile radio with L g = L h = 4. The small step size a is 
chosen to be 0.002. The performance of linear prediction relies on the constraint 
value y due to its impact on the MSE behavior. Throughout this thesis, y is chosen to 
be 1.5. Figures 4.18 and 4.19 illustrate the behavior of this algorithm in a white noise 
situation. For simulation, use equations (3.31)— (3.34). 
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Figure 4.18 Eye diagram of equalizer output, using LP method with the raised-cosine 

channel. 
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Figure 4.19 Mean-square symbol error of LP method. 

From Figure 4.18 and 4.19, we see that this algorithm works well to equalize 
the received noisy signal. The MSE converges to 0.06786 at index sample n = 
50000. The problem with this algorithm is it has slow convergence. The amount of 
rotation and magnification in the eye diagram (see Figure 4.18) is a function of 
/i (0) = (0.52 —jO.72) which is the leading coefficient in h . The angle of rotation 
is 54.16°, and the magnification is |/i (0)| = 0.88. 
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4.2.7 Experiment 4b 

The second class of channel model is simulated in this experiment with using 
the same method as in experiment 4a. Our assumptions are still the same. Figure 4.20 
shows the equalized eye diagram, while Figure 4.21 depicts the mean-square symbol 
error across 50000 symbols. 
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Figure 4.20 Eye diagram of equalizer output, using LP algorithm with the second 

class of channel model. 
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Figure 4.21 Mean- square symbol error of LP algorithm. 

The algorithm works very well as we see from Figure 4.21 and the MSE 
converges to 0.0608 at time sample (50000). The performance is slightly better than 
that in experiment 4a. However, we are still facing the slow convergence problem. 
The leading parameter /i (0) = (0.5498 — jQ. 3174). The angle of rotation is 30°, 
and the magnification is |/i (0)| = 0.635. 
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4.2.8 Experiment 5a 

In this section, we use a symbol sequence that is assumed to be uncorrelated 
and generated from a 16- QAM constellation with -1.5, -0.5, 0.5 and 1.5 levels. The 
RELS algorithm, which is first presented in [22], is simulated using the following 
(third class) IIR channel model 

A{q~ 1 ) = l-Q.Qq- 1 + QAlq~ 2 . (4.9) 

Fiq- 1 ) = 1 + (0.5 - 0.60<7 -1 . (4.10) 

Biq' 1 ) = 0.5498 - 0.3174i + (0.4342 - 0.45960q _1 + (0.0131 - 

0.0919i)q" 2 + (-0.0275 + 0.1592i)q" 3 - (4.11) 

Ciq' 1 ) = 0.0414 - 0.0717i + (0.7168 - 0.3407i)q~ 1 + (0.0745 - 

0.3636i)q" 2 + (0.0685 - 0.1861i)q" 3 . (4.12) 

Our assumptions with respect to additive noise are still the same. Figure 4.22 
shows the received symbols x 1 (i), while Figure 4.23 shows the equalized symbols 
eye diagram. As performance measures, the mean-square symbol error for RELS 
algorithm across 3000 symbols is estimated as shown in Figure 4.24 using the 
following equation 

mse(n) = ^£?=ol*i(i + 1) - 9(i + 1) - b w(i + 1)| 2 . (4.13) 

For simulation, use equation (3.49)-(3.53). 
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Figure 4.22 Eye diagram of received signal for 16 QAM modulation. 
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Figure 4.23 Eye diagram of equalizer output, using RELS algorithm with the third 

class of channel model (IIR). 
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Figure 4.24 Mean square symbol error of RELS algorithm. 



4.2.9 Experiment 5b 

In this section, we run the algorithm with the presence of the modeling error 
in the channel dynamics, and we use the same polynomials as in Experiment 5a, 



B(q-i)F(q-i) 



*i(0 = — A r _!-v w(i) + v x (i) + 8 x (i), 



Ciq-^Fiq- 1 ) 



x 2 (0 = A( _!•> w(i) + v 2 (i) + 6 2 (j0. 



(4.14) 
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where 1^(0 and v 2 {i) are white noise with zero mean and variance=0.14. It is 
important to mention that u k (i),k = 1,2, and w(i) are independent sequences. Let 
8 k {i), k = 1,2, are modeling errors given by 



*i(0 



1- 0.5cri W(0 aUd ^ (0 = l + 0.7^ w(a (4 " 15) 



where 5 is a parameter defining the size of the channel modeling error. Figures 4.25 
and 4.26 illustrate the scatter plot of the equalized symbols and sample mean-square 
error, respectively, with setting 8 — 0.1 in (4.15). 
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Figure 4.25 Scatter plot of equalized symbols, using RELS algorithm with the 
presence of modeling error using (IIR) channel. 
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Figure 4.26 Sample mean square error of RELS algorithm. 

This algorithm shows a very good performance of equalizing the received 
signal with and without presence of noise and modeling error as illustrated in both 
Figures 4.24 and 4.26. The MSE converges to 0.04598 at index sample n = 3000 in 
experiment 5a, whereas the MSE of experiment 5b converges to 0.06469 at the same 
index sample. Therefore, it is clear that the algorithm provides a degree of robustness 
with respect to receiver noise and modeling error [22]. 

The comparison between all of these algorithms is done with using the second 
class of channel model. Figure 4.27 demonstrates that the performance of RELS 
algorithm is the best since its MSE has the lowest value at the time sample (n = 
3000). 
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Figure 4.27 Mean-square symbol error, Comparison for those five algorithms with the 

second class of channel model. 
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CHAPTER V 

CONCLUSION 

This thesis has proved how adaptive blind algorithms (RLS, cyclic LMS, 
LMS, LP and RELS) can equalize the wireless communication channels through 
removing inter-symbol-interference (ISI). These techniques have proven to be quite 
effective and powerful to combat ISI effect. ISI is mainly generated in dispersive 
channels such as Radio and Mobile wireless channels. In a particular case, the mobile 
cellular communication (multipath propagation of the transmitted signal) suffers from 
severe ISI. Therefore, an adaptive filter as an equalizer can be placed at the receiver 
to compensate or equalize the dispersion occurred during transmission. This filter or 
the equalizer is a device positioned at the receiver to alleviate the effect of ISI and 
thus the transmitted symbol sequence can be recovered. 

Channel equalization methods used in this thesis rely only on the statistical 
behavior of the received signals (Second order Statistics) in order to estimate the 
transmitted sequence without requiring the knowledge of the channel characteristics. 
The transmission environments used were a raised-cosine SIMO channel model 
realized by fractionally sampled (FS) FIR filter, SIMO channel model with constant 
coefficients over period realized by (FS) FIR filter and SIMO channel model realized 
by (FS) IIR filter. The latest channel model was performed with and without 
modeling error. 
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It has been shown that the algorithms used in this thesis perform very well, 
and they are very efficient and robust with respect to the channel distortions through 
examining their mean square symbol error (MSE) performance via computer 
simulations. Also, it is shown that the RELS algorithm has an advantage over the rest 
of algorithms since it has the fastest MSE convergence. 

For future research, our system might be extended to MIMO channel model 
case, and we should study and examine our techniques performance. Also, we may 
investigate the capability of adaptive blind channel equalization algorithms with using 
time-varying transmission channel and see which one provides computationally 
efficient implementations with and without presence of noise. 
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